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In The Specification: 

Replace the following paragraphs as follows: 

[0018] In this example two images A1 and A2 of the same object or scene are to be 

imaged matched A2 to A1 . For every pixel of A1 and A2 % the following relation holds: A1 ■ A2 * 
(A1/A2). By difforontiation of th e logarithm of above e quation, th e contract function G(.)a\ a 
givon l ocation ic donotod by : C(A1 ) - G(A3) i C(A1/A£). Ac wil l bo further docor i bod below, tho 
imag e division A1/A2 may optionally bo rogu l ar i aod relative to tho i mago to bo matohed A1 i n 
ctop 52 whon tho imago qua l ity i o not good, o.g., noisy. Various typoo if rogu lafi zat i on may b e 
porformodj Rogular i zat i on wil l b e furthor doscribod bolow. 
A1 = A2*(A1/A2V 

{0 Q . 1 8 ] I n ordor to catiofy C(A1) - G(A2) i n tho abovo e quation, G(A1/A2) » Q. A w ell known 
way to dooroas e tho contrast is to l ow pace filtor tho ratio A1/A3 ae shown in ctop 61. Thoroforo 
in-a t e p 56y contrast matching output equ a tion for tho two imagos A1 and A2 i c thuc: AW-^Ag-^ 
LPF(A1/A2), whore At^ i o tho contrast matchod vorcion of A2 with rospeot to A1 and LPFQ i s 
a l ow pass fi l tor funct i on. Th e l ow paoo f il t e r function io further d oc or i b o d bo l ow, 

[00191 Bv differentiation of the logarithm of above equation, the contrast function CD at 

a given location is denoted by : 

C(A1) = C(A2UC(A1/A2). 

[0020] For mu l tiplo (N) imag es , — Lot A1, A2, AK~, AN bo tho N imag es und e r 

consid e ration (K r <AQ and oaoh of thoso imag es ar e to bo matched to th e e amo roforonoo i mag e 
A1. By e xt e nd i ng tho abovo l og i c to any of N i magoo, cay imago K r tho g e neral relationship 
e xis ts, - A W = AK * LPF(A1/AK) whoro Al^ - to - th e contract matchod vorc i on of 4* with r es p eet 
to A1, Thuc, a gonorafe e d contra s t match i ng has boon achiev e d si nc e , C(A1) = Cf/l W m ■ 
£{Al M)-...aC(A1 M ^ 



f 00201 As will be further described below, the image division A1/A2 may optionally be 

regularized relative to the image to be matched A1 in step 52 when the image quality is not 
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good, e.g.. noisy. Various types if reoularization may be performed. R eoularization wit* be 
further described below. 

[00211 In order to satisfy CfAD = C(A2) in the above equation. C(A1/A2) - 0. A well 

known way to decrease the contrast is to low pass filter the ratio A1/A2 as s hown in step 54. 
Therefore in step 56. contrast matching output equation for the two images A1 an d A2 is thus: 
A1m» =A2«LPF(A1/A2). 

where Al^g is the contrast matched version of A2 with respect to A1 and LP Ff.) is a low pass 
filter function. The low pass filter function is further described below. 

[0022] T o - summarize an i mago A2 has to b e matchod to another imago A 1 of the oarno 

scon e /obj e cts to obtain tho matchod imago A1x& uoing tho ro l ation: Al^ a _ » A2 » LPF(A1/A2) 
whoro LPF i o a l ow pass fi l t e r function. Proforably tho l ow pass filt e r funct i on i o a boxcar filter 
and tho param e t e rs of tho f il t ef- aro application spec i f i c. — For gonora l app l ications, th e f i lt e r 
kornol longth is on e tonth tho longth of tho imago (assuming a cquaro i mago and square 
kornol). Furthermore, in praotioo, the abovo oquation may nood to bo mod i f ie d i n ordor to avo i d 
nolco amp li fication during i mago d i vis i on. Rogular i aation may bo porformod i n a numbor of 
me thod s to provont nois e amplif i cation during i mago diviG i on. Th e i mago division ratio ha s a 
numorator A1 and a donom i nator A2. On e mothod to r e gular i z e imago divis i on io to add a emai l 
oon s tant to tho denominator, i. e ., denominator b e com e s (A2 i c). wh e r e a s an oxamp l o, c » 1 .0. 
Thus th e e quation booomos » A2 * LPF(A1/(AS > o)). 

f00221 For multiple (N) images. Let A1. A2. .. AK... AN be the N images under 

consideration (K<N) and each of these images are to be matched to the same reference image 
A1 . By extending the above logic to any of N images, sav image K. the general relationship 
exists. 

A1 ^ = AK*LPF(A1/AK) 

where A1^ is the contrast matched version of AKwth respect to A1. Thus, a generalized 
contrast matching has been achieved since. CfAD = C(A1^) - ... = CfAl^) = ...= CfAluy l 
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{0024^100231 The choice of parameters in the low pass filter function essentially determines 
the scale of contrast matching obtained. Various types of low pass filters may be used. For 
example, a boxcar filter with a single parameter may be used. A boxcar filter smoothes an 
image by the average of a given neighborhood of pixels. It is separable and efficient methods 
exist for its computation. Each point in the image requires just four arithmetic operations, 
irrespective of the kernel size. The length of the separable kernel is variable and depends on 
the scale of contrast matching desired. For example, if the kernel size is about one tenth of the 
image size, assuming a square image and a square kernel, excellent global contrast matching 
of images is obtained. On the other hand, using too small a kernel size produces undesirable 
blobby patterns in the matched images. Therefore, a reasonably large kernel should be used to 
avoid any perceptible artifacts using this method. 

[0024] Anothor me th od for rogularization io to roplaoe tho ratio (A1/A2) by a regular i z e d 

ratio givon by (A1 * A2/(A3 * A2 i S)) t wh o ro ac an oxampl e , S « 1.0. Thuo tho oquation 
bocomoo AW y A2 * LPRA1 * AQ/(A2 * A2 I 5)). 

[0024] To summarize, an image A2 has to be matched to another image A 1 of the same 

scene/obiects to obtain the matched image Ah *» using the relation: 
A1^ = A2»LPF(A1/A2) 

where LPF is a low pass filter function. Preferably the low pass filter function i s a boxcar filter 
and the parameters of the filter are application specific. For general applica tions, the filter 
kernel length is one-tenth the length of the image (assuming a square ima ae and square 
kernel!. Furthermore, in practice, the above equation mav need to be modified i n order to avoid 
noise amplification during image division. Popularization mav be performed in a number_of 
methods to prevent noise amplification during image division. The image division ratio has a 
numerator A1 and a denominator A2. One method to regularize image division is to add a small 
constant to the denominator, i.e.. denominator becomes (A2 + 1). where as an e xample, e - 1.0, 
Thus the eouation becomes 

A1^ = A2 * LPRA1/(A2+£tt. 

[0023] r00251 Of course, if no regularization is to be performed, £ would be 0. 
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[00261 Another method for requiarization is to replace the ratio (A1/A2) by a regularized 

ratio given by (A1 * A2/(A2 * A2 + S)\ where as an exa mple. 5 = 1.0. Thus the equation 
becomes 

A1 », = A2 « LPRA1 * A2/(A2 - A2 + 5)). 

teoaa r00271 When a number of images A2, AK,... f AN have to be matched to a single 
image A1, the above process may be performed in a pair wise fashion to obtain A7m2>..., 
A1mK"uA1mm 

{08881 [00281 While the invention has been described in connection with one or more 
embodiments, It should be understood that the invention is not limited to those embodiments. 
On the contrary, the invention is intended to cover ail alternatives, modifications, and 
equivalents, as may be included within the spirit and scope of the appended claims. 
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